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f open resonators to measure electro-optical 
iielectric samples and surface resistivities of 
films has been recently extended into the 
wavelength region. A primary study was devised 
with a scanning hemispherical open resonator 
•=ze at these frequencies. In their work, a 
- v is fed by a far-infrared laser ope rat ing at 
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Resonator methods make use of the chance in cavity 
quality Q (discussed later) and change in resonant length to 
measure metallic surface resistivities and complex dielectric 
permittivities. Typically, at wavelengths larger than 10 mm 
(< 30 GHz) , classical closed cavities have been used in these 
measurements. At some short-millimeter wavelengths, tne small 
size recruited for these microwave resonant closed cavities has 
limited the ability to mechanically construct cavities. As we 
approach still shcrter-millimeter wavelengths, it first 
becomes inconvenient and* eventually impractica. to prepare 
samples for closed cavities. An open hemispherical resonator 
has been proposed (see Ref. 1) as an alternative to tne 
waveguide in the investigation of the characterization or 
material properties at higher frequencies, in particular, m 
the region between 337 GHz and 2.7 THz . Since these 
frequencies lie between the optical and microwave regions, 
the open resonator is an attempt to combine useful techniques 
from both regimes . 

In the aforementioned studies, the measurements of the 
surface resistivity or dielectric constant were made by 
placing samples on the planar mirror in the hemispherical 
resonator. The energy dissipated or lost in tne cavity is 
mainly in the conducting surfaces 3 while total stereo energy 

_ 4 

defends on the Geometrical dimension of the resonator, tones 
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adapted the Gaussian bean formula for use with a hemispherical 
resonator and described theoretical calculations of- cavity Q, 
which can be obtained from the usual electromagnetic fields 
theory . 

In order to theoretically predict accurate Q values, it 
is important to be able to calculate the exact total energy 
stored in an open resonator. In the hole-coupling approach in 
the experiments performed by Cook; et ai. (see Ret. 1) wnere a 
laser feeds the scanning resonant cavity through a small hole 
centered in the spherical, mirror, the energy loss through a 
coupling hole should also be considered in measuring the Q 
factor of the cavity. They estimate a 2% error in calculating 
total energy in their work. 

The purpose of this work is no extend resonator theory 
into the region in which the planar mirror is quite smal_. 
P.esults of the theoretical description are then enter. ced to 
resonator design and experimental arrangements as discussed in 
further sections of this work. Finally, a discussion os 
dielectric measurements for small samples is included as a 
specific application of this work. 



1 . Review of open resonator theory 


Before discussing specific applications, it is first 
necessary to review existing theories of open resonators. An 
open-resonator system uses two metal reflectors to form a 
resonant structure similar to that of a Fabry-Perot etalon. 
In making measurements such as surface resistivity or 
dielectric constant, the open resonator is of the 
hemispherical type and consists of one spherical and one 
planar mirror separated by a distance slightly less than the 
radius of curvature of the curved mirror as shown in Figure 1. 
This figure shows a typical open resonator with dimensions 
scaled for the submillimeter-wavelength region. Figure 1 also 
shows a Gaussian field scored in the cavity. This field is 
characterized by a beam waist w D . 



The hemispherical rype of open resontor 
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This structure has several advantages. First, it is 
easily aligned. In fact, small misalignments can be ignored. 
Second, the experimental setup is easily arranged. Also, the 
beam waist is on the planar mirror, which serves as a sample. 
Furthermore, this system has a high cavity quality factor, or 
Q, defined as 2 n times the ratio of the time-averaged energy- 
stored in the cavity to the energy loss per cycle. The mode 
configurations of open resonators are characterized by the 
Kirchhoff -Fresnel diffraction theory 5 of electromagnet ic 
radiation in which the electric field (or the magnetic field) 
is represented by a single scalar runction E as describee 
below. In view of the importance of this theory, and its 
relevance to this work, the derivation will be sketched here. 
The following derivation is based on one descripLJ.ons “a- ven 
Kogeinik and Li 5 and Maitland and Dun.r. 

1.1 Spatial field distribution in open resonators 


The scalar theory 
(Helmholtz) : 


leads to the scalar wave equati< 


v i E + k‘ ; E = 0 (1.1) 

where k=2 n/X is the propagation constant in free space ana 
E=E(x,y,z) is the scalar field, describing a wave travelli-.- 
in the + z direction. (For convenience time variation has bee.. 
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icrnored.) A travelling wave solution to the electromagnetic 
wave equation should be of the form: 


E (x, y, z) =V F (x, y, z) exp (- jkz) 


(1.2) 


where is an unknown amplitude function to be deter- 
mined by substituting Eq. (1.2) into Eq. (1.1) . The complete 
expression for travelling solutions are found (see Ref. 5) to 
take the form : 


E (x, y, z) =E a (w 0 /w) [H,, Cv2x/w) H n (V2y/w) exp (-r 2 /w) ] 
• exp (- j (kz+<J)) - jkr 2 /2R] 


The terms in the square brackets in Eq. (1.3) describe 
the transverse amplitude variation ot tne beam.. H n a. c 
Hermite polynomials which describe variations or tne mode 
in the x and y directions, respectively, r is the transverse 
distance from the z axis. 
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Figure 2 Geometrically shows trie relationship between w 
and w Q in Eq. (1.3). w=w(z) is the radius parameter of the 
beam given by 

w 2 =w 2 [ 1+ (2z/kw 2 ) 2 ] (1.4) 

where w Q is the minimum radius at the waist of the beam and is 
determined by boundary conditions. The distance of the beam 
from the origin (z = 0) is -z which can be expressed in terms of 
the parameters of the cavity 5 . For z x and z 2 in Figure 2, we 
consider the Gaussian beam as the primary mode or a resonant 
cavity formed from two spherical mirrors (radii of curvature 
r x and R 2 ) with a distance D apart. They are given by 

Zi=- (R 2 -D>D/ (R : +R 2 -2D) (1.5) 

z 2 = (R,-D) D/ (R,+R 2 -2D) (1-6) 

The term <j) = (})(z) describes a pnase change on the axis of the 
beam given by 

d)= (1+m+n) tan' 1 (kw 2 /2z) (1.7) 

In the argument of the final exponential or Eq. (1.3) , R— R(z) 
is the radius of curvature of the wavefront. It is given by 


R=z [1+ (kw 2 /2z) 2 ] 


( 1 . 8 ) 
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R becomes infinite (i.e. describes a plane wavefront) at the 
waist. At the reflectors, R is equal to the radius of 
curvature of the reflectors R x and R 2 . This condition 
determines w Q , 


(kw^) 2 =4 (R : -D) (R 2 -D) (R!+R 2 -D) D/ (R! + R 2 -2D) 2 (1.9) 

In the hemispherical cavity, R 2 is slightly larger than D and 
R x goes infinite. Eq. (1:9) reduces to 

(kw 2 ) 2 =4 (R 2 -D) D (1.10) 

Mathematical superposition of forward and reverse 
traveling beams of similar mode and the fitting of tne 
resulting standing wave to the boundary conditions at the 


reflectors 

gives the field 

description of 

: the 

TEM mode 

resonances , 

where m,n=0,l,2., 

. . refer to vari 

ation 

in the x, y- 

directions 

and specific transverse modes. 

q is 

an integer 


typically much greater than 1 and describes longitudinal 
modes . 

For the cavity to resonate, the phase shift when the beam 
travels from one mirror to the other must be an integral 
multiple of ft because only in this case will a field be 


established inside 

the cavity with a 

well 

def i 

nee. 

phase 

structure. Using t 

he phase term of Eq. 

(1.3) 

for 

the 

phase 

shift per transit. 

along the z-axis 

( r=Q ) f 

we 

have the 
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following condition for resonance: 

kz 2 + (t> (z 2 ) - [kz t +<|> (Zi) ] =q7t (1.11) 

where q is an integer which equal to the number of the 
standing waves in the cavity. Inserting Eq. (1.5), (1.6), and 

(1.7) into Eq. (1.11), one obtains 

kD=qm+ ( 1+m+n) tan'V D (R\+R 2 -D) / (R^-D) (R 2 -D) (1.12) 

The resonant frequency of the resonator for the TElL„ q mode is 
given by 

f^ q =c [qm+ ( 1+m+n ) tan' 1 ’/ D (R,+R 2 -D) / (R : -D) (R 2 -D) ] / 2 tcD ( i . i o) 

and the resonant length D is 

D=c [q7C+ (1+m+n) tan -1 -/ D (R 2 +R 2 -D) / (R 1 -D) (R 2 -D) ] /27Cf mnq (1.14) 

Since the radius R : goes infinite in the hemispherical 
resonator, from Eq. (1.13) one therefore obtains the 
resonant frequency: 


f anq 


-c [ <q7t + ( 1+m+n) can'V D/ (R 2 -D) ] /2kD 


(1.15) 
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1.2 Resonator without sample 

By definition the cavity quality Q can also be written as 

Q=f (energy stored) / (power dissipated per radian) (1.16) 

where f is the radiation frequency. The total electromagnetic 
energy density is denoted by 

u= (£ 0 E 2 +^ 0 H 2 ) /2 (1.17) 

where e 0 and |i 0 are electric permittivity and magnetic 
permeability in the vacuum, E is the electric field intensity, 
and H is the magnetic field intensity. For any transverse 
electromagnetic or TEM wave, the electric and magnetic energy 
densities are equal, 

E a E 2 / 2=ll 0 H z / 2 (1-18) 

and the total energy density is therefore £ C E 2 . It losses due 
to diffraction or coupling hole are neglected, then the energy 
stored in an empty resonator is expressed as 

W c =eJ v <E 2 >dv (1.19) 

where dv in cylindrical coordinates is expressed as 27irdrct, 
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V is the total volume of the cavity, and the time average term 
<S 2 > is E'E’/2. Since only the fundamental mode is considered, 
m=n=0 and both Hermite polynomials reduce to unity. The term 
exp (-r 2 /w 2 ) in Eq. (1.3) describes the dominant Gaussian 

profile of the beam. Thus 

<E 2 >=[E 2 w 2 (exp (-2r 2 /w 2 ) /w 2 (z) ) ] /2 (1.20) 

Substituting Eq. (1.20) into Eq. (1.19) gives 

W Q = [e o E 2 w 2 J v (exp (-2r 2 /w 2 ) /w 2 (z) ) dv] /2 
=t£ 0 E 2 w 2 J° [J- (exp (-2r 2 /w 2 ) /w 2 (z) ) rdr] d 

In Eq. (1.22), the limits of integral of r go from zero to 
infinity, i.e. the cavity is considered to have infinite 
dimensions. Setting x=r 2 /w(z) 2 , we have 

dx=2rdr/w (z) 2 (1.23) 

and W 0 =7t£ 0 E 2 w 2 j° [j“exp (-2x) dx/2] dz (1.24) 

= [7t£ 0 E 2 w 2 |°dz] /4 (1.25) 

If the separation of the mirrors in the resonator is denoted 
by D, then 


( 1 . 21 ) 

( 1 . 22 ) 


W 0 =7Tw 2 £ o E 2 D/4 


(1.26) 
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Eq. (1.26) agrees with the expression given by Siegman 7 and 
differs from Jones (see Ref. 4) by a factor of 2. 

Now let us consider the energy losses in the metallic 
mirrors. The ratio of displacement current density dD/dt to 
conduction current density J f in a conductor can be estimated 
as 8 

I (8D/9t) / Jf I =27t£ 0 f/C (1.27) 

We may set cr=10 7 ohms/meter for a good conductor, then 

| (3D/3t) /J f | =10' 17 f (1.28) 

where f is the radiation frequency. The displacement current 
in a good conductor is therefore negligible at any frequency 
lower than 10 15 hertz. With this approximation, the Maxwell 
curl equation becomes 

E c =VxH c /a (1.2 9) 

where vectors E c and H c represent electromagnetic waves in the 
conductor. If n is the unit normal outward from the conductor 
and £ is the normal coordinate inward into the conductor, 
then the gradient operator in Eq. (1.29) can be written as 


V=-nd/d£ s 


(1 . 30) 
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Since inside the conductor H is parallel to the surface, the 
solution for H c is 9 

H c =H 0 exp (-£/8) exp (i^/5) (1.31) 

where H 0 is the tangential magnetic field outside the surface 
and 8 is the skin depth defined by 

5= (l/7tf aji 0 ) :/2 (1.32) 

Inserting Eq. (1.30) and Eq. (1.31) into Eq. (1.29), we 
obtain 


= (-nx3H c /9£) / a 


(1.33) 

={ (l~i) exp [ (i-1) 

i ^/5] nxH 0 } /g5 

(1.34) 

=R (1-i) exp [ (i-1) 

\ ^/8] nxH 0 

(1.35) 


where R=l/cr5 is defined as the surface resistivity of the 
conductor. The time average power absorbed per unit volume is 

p= ( j f - E c *) /2 (1.36) 

=OE e - E c 72 (1.37) 

=aR 2 exp (-2^/5) H 2 (1.38) 

where indicates the complex conjugate. The total power 


loss is 
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?=/ L pAd^ (1.39) 

where A is the area of the beam cross section rcw 2 ana L is the 
thickness of the conductor. Because L is much larger than 5, 
then Eq. (1.33) becomes 


?=J“pJtw 2 d^ (1.40) 

=RH 2 7tw 2 / 2 (1.41) 

=RmwfEo/2 Zo‘ (1.42) 

where Z 0 = (|l o /e 0 ) 54 is the impedance of free space. The net 
power P t loss in the resonator is given by 

P c =P c +P p =7tw 2 E 2 (R c +Rp) /2Z 2 (1-43) 

=TCw 2 E 2 R c (M+l) /2 Zc (1.44) 


where M is Rp/R*. and subscripts p and c indicate the planar 
mirror and concave mirror, respectively. From Ec. (1.16), the 
Q factor of the empty resonator is finally given by 

Q 0 =f W 0 /P t (1.45) 


= fD£ 0 Z 2 /2R c (M+l) 


(1.46) 



2 . Small samples consideration 


We now extend the theory to cover situations where the 
planar mirror is quite small. When making measurements of 
surface resistivities of metallic materials, samples are used 
as planar reflectors to form one mirror of the empty resona- 
tor. Some samples are small due to expense and availability 
so that modifications of resonator theory are required. The 
results can then be applied to a resonator loaded with small 
dielectric samples. The complex dielectric constants of these 
samples can be extracted from such measurements. 

2.1 Spurious modes in the hemispherical^ cavity 

The theory used by Jones 10 applies only to the Gaussian 
TEM ooq modes supported in the cavity. Resonant frequencies 
other than the TEM ooq modes were given by Eq. (1.15) . .Those 
resonances were attributed to spurious (higher-order) modes 
and had been ignored since they were sufficiently separated 
from the Gaussian mode. However, the experiment performed by 
Cook et si. (see Ref. 2) has shown that the spurious modes are 
not negligible as the sample radius is enlarged. These 
spurious modes may lead to confusion in the resonator cavity, 
where higher-order modes from frequencies other than the 
Gaussian resonant frequency could be supported. To prevent 
the introduction of spurious modes within the resonant cavity, 
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it was suggested (see Ref. 2) that the radius or effective 
sample size of the planar mirror could be reduced into a 
working radius slightly larger than the Gaussian beam spot 
w Q . This size, while large enough to support the desired 
Gaussian TEM ooq mode, is small enough to prevent the occurrence 
of higher-order modes. However, it is not clear that existing 
resonator theory is adequate for this application. We next 
address this question. 

2.2 Modification of the theory 

To proceed further we must investigate two situations: 
(1) The small sample does not affect the electric field 
distribution in the cavity so that the field of the fundamen- 
tal modes is an exact Gaussian beam, and (2) the consequent 
possibility of diffraction losses can still be neglected. 

In an open resonator formed by a hemispherical cavity, 
fundamental TEM ooq modes have Gaussian profiles and the 
electric field varies in the radial direction as 

E (x, y, z) =E 0 (w 0 /w ) exp (-r 2 /w 2 ) exp (- j (kz+<J)) - jkr 2 /2R) (2.1) 

Since the small sample serves as a planar mirror, assumption 
(1) is valid as long as the sample is mounted without mis- 
alignments in the cavity and large enough to support the 
Gaussian Mode. The beam radius w given by Eq. (1.4) varies 
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parabolically along the axis having a minimum value, w Q , at 
the waist (z=0) of the resonance. In experimental arrange- 
ments, the sample is placed in the position of the beam waist 
in the hemispherical cavity and its size is sufficiently 
larger than the beam waist w Q so as to prevent "leakage". 
Assumption (2) is also true since the Q factor of the resona- 
tor due to the diffraction loss was estimated to be over 10 7 
in these experiments 1 * while the measured Q factors (given in 
Table 1 of section 3.5) have values on the order of 10 5 . The 
diffraction loss in the cavity is therefore negligible 
compared to the conductor losses. 

Reported measurements from Jones (see Ref. 10) are based 
on the hemispherical cavity in which the sample diameters were 
often in the range of 70-80 mm. These diameters are much 
larger than those of the Gaussian beam in the resonator. The 
Jones' theory assumes that the electric field is negligible 
for transverse distances much larger than the beam waist w 0 . 
When the sample is of a small radius compatible with the beam 
waist w 0 , this assumption no longer holds. 

To accurately calculate the Q factor of the cavity, we 
first modify the stored energy expression (1.22). Consider 
the hemispherical cavity shown in Figure 1. The total energy 
stored in an open infinite cavity is given by Eq. (1.19) . For 
the empty resonator with a small planar mirror, the integral 
of r does not go to infinity in Eq. (1.22) . Thus 
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W 0 =£ 0 E 2 w 2 J° [J; (exp (~2r 2 /w 2 ) /w 2 (z) ) rdr J dz (2.2) 

= [7I£ 0 E 2 w 2 J° (1-exp (-2r 2 /w 2 ) ) dz ] /4 (2.3) 

When r is much larger than w, we find W Q equals r£ 0 E 2 w 2 D/4 as 
shown in Eq. (1.26). In the case of small samples, the 
transverse distance r in the resonator is just a little larger 
than the beam radius w. By setting r=bw (b>0) , where b is a 
scaling constant, from Eq. (2.3) we can see that 

w'=tc£ 0 E 2 w|D [1-exp (-2b 2 ) ] / 4 
or W„=W = [ 1-exp (-2b 2 ) ] 

Eq. (1.46) is then changed into 

Qo=Q 0 [1-exp (-2b 2 ) ] (2.6) 

The theoretical curve from Eq. (2.5) is shown in Figure 3. In 
this figure, we see that the ratio of the corrected energy W 3 
to the uncorrected energy w o is very close to 1 as long as r 
is not smaller than 1.5w. 

The importance of this result can be summarized as 
follows; 

(1) At the point of r/w=1.5 in Figure 3, the corrected energy 
W c is 0.9889W,. Thus, for the sample size slightly larger than 
1.5w,, the energy leakage can be ignored while the spurious 
modes have been effectively repressed. 


(2.4) 

(2.5) 
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as a function of planar mirror sizes 
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(2) This has the practical advantage of permitting samples of 
smaller radius to be used in measurements since the beam 
radius w has its minimum value w Q at the planar mirror in the 
hemispherical resonator. 

(3) The Q-factor of cavities containing small samples can be 
more precisely predicted since we know the exact energy stored 
in the cavity. 

Before one can optimize the design of an open resonator 
to suit it to a particuiar application, it is important to 
note the experiment result given by Cook et al. (see Ref. 2) . 

The experiment is performed at 337 GHz with R 2 =135.9 mm and 
D=135.5 mm. Their work shows that spurious modes are prevent- 
ed when the radius of the sample is smaller than 5 mm. The 
beam waist can be obtained by Eq. (1.10) . It gives w 0 =1.5 mm. 
Using the theoretical results for the consideration of no 
energy leakage in the cavity, we can select the sample radius 
in the range of 1 . 5w a (=2.3 mm) to 5 mm. 

It should be pointed out that Eq. (2.6) must be used to 
qive precise measurements when the sample radius is smaller 
than 1.5w a . However, the diffraction loss in such small 
samples may not be neglected. The magnitude of the diffrac- 
tion loss can be estimated (see Ref. 11) in determining the Q 
factor of the resonator. This diffraction limits the effec- 
tive size of planar samples to about 1.5w 0 . 
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2.3 Resonator containing dielectric sample 

The open resonator used in making precise measurements of 
permittivity and loss angle (discussed later) for a range of 
dielectric materials is mainly the hemispherical type devel- 
oped by Jones (see Ref. 10) . A "loaded" hemispherical 
resonator consisted of a planar mirror and a concave mirror is 
shown in Figure 4. 



Fig. 4. Hemispherical resonator 

A=concave mirror 
Biphase front 

Oplane mirror at resonant posicion with sample 
I=plane mirror at resonant posicion without sample 
F = s amp 1 e 
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After making the empty cavity into resonance, a dielec- 
tric sample of thickness t, with refractive index n, is loaded 
onto the planar mirror as indicated in Figure 4, a change p in 
length of the resonator is required to bring the system back 
into resonance. Following the theory describe by Jones (see 
Ref. 4), the stored energy expression (1.26) is 

W.-TCwifioE* (tA+d) / 4 

where A=n 2 / [n 2 cos 2 nkt-<J) c ) +sin 2 (nkt-dn) ] 

and (|) t =tan' 1 (2t/nkw 2 ) 

Similarly, the Q-factor from Eq. (1.46) becomes 

Q L =fE 0 Z| (tA+d) /2R C ( M A - 1 ) (2.10) 

Considering the resonator with the small sample and following 
the calculations used for an empty resonator, Eq. (2.5) can be 
written as 

W^=W L [ 1 -exp ( -2b 2 ) ] (2.11) 


(2.7) 
( 2 . 8 ) 
(2 . S) 


Correspondingly, Eq. (2.6) becomes 


Q^=Q; [1-exp (-2b 2 ) ] 


ratio of Q[ to Q~ is given by 


(2 . 12) 


The 
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Ql/Qc=Q-,/Qo=(M+1) (tA+d) / (MA+1)D (2.13) 

Ql = C; O'-l) (tA+d) / (MA+I)D (2.14) 


This last equation will be used for calculating the loss angle 
later . 



3. Using an open resonator for empirical dielectric 
measurements 

Determination of dielectric properties of materials is a 
natural application for open resonators since easily measured 
parameters of the resonator, the resonant irequency f, 
resonant length D, and the Q factor, are simply related to tne 
permittivity and loss of the included media. This chapter 
will cover two applications of the open hemispherical resona- 
tor . 

3.1 Electromagnetic quantities 

The complex refractive index n is derived from complex 
dielectric permittivity £ of Maxwell 7 s equations so that 


The real 


and imaginary parts of n are, 


by definition 


n=n-iK=n-i (ac/4rcf) 


(3.2) 


where k is 
coefficient 
of light in 
complex diel 


the absorption index ac/47tf, a the absorption 
in cm -1 , f the frequency in hertz, c the velocity 
vacuum, and n the real retractive index. ir.e 
ectric permittivity £ has real and imaginary parts 
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expressed as 


£=£ ' -i£" 


(3.3) 


Then, our definitions provide us with the simple relationships 
between the fundamental optical quantities, a and n, and the 
dielectric quantities £' and £", as follows 

e' =n 2 -K : =n 2 - (acV4jrf ) 2 (3.4) 
e"=2nK=nac/27tf (3.5) 

The term loss tangent, or as is commonly expressed, tcnS, is 
the ratio of the imaginary part £" to the real part £' of the 
dielectric permittivity 

tan5=£" /£' (3-5) 


Note that the expression for the dielectric constant E ' is not 
a constant. It contains a term inversely dependent on 
frequency. It is net a constant, that is, unless the second 
term that contains the frequency can be neglected. This is 
usually the case at optical frequencies where the denominator 
of the second term is very large. At millimeter wave frequen 
cies, however, we cannot drop the second term unless the 
absorption coefficierc ct is sufficiently small c^rn^a-sd wiv.:, 
the refraction index n sc that the real part of tne dielectric 
permittivity £ ' is effectively independent cf frequency. 
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Specially, if the absorption coefficient is less than unity at 
300 GHz, then £' is truly a dielectric constant and the 
material is "low loss" in chis case. 

3.2 Open resonator method 

Resonator methods make use of the change in Q and in 
resonant length to provide £' and £". A convenient form of 
resonant structure is the hemispherical type that employs one 
concave and one planar mirror as shown in Figure 4. This 
geometrical consideration may allow measurements to be made on 
small diameter samples because the beam waist spot is on the 
planar mirror in the cavity. The sample placed on the planar 
mirror can be a liquid or a flat solid. Cullen and Yu : ", 
considering a sample of thickness t, have applied the beam 
wave theory to the open resonator containing a dielectric 
sample and have derived equations from which the refractive 
index n and tangent of loss angle 8 can be obtained. These 
equations are 


tan (nkt-<b c ) /n=-t an (ka-<j> d ) (3.7) 

d=D-t-p (3.8) 

(ji^tan* 1 (2t/nkvq) (3.9) 

(j) d =tan' 1 [2 (d+t /n 2 ) /kvil - tan -1 (2t / n 2 kw|) ( 3 . 10) 

kw 2 =2 [ (d+t/n 2 ) (Rg-d-t/n 2 ) ] 1/2 (3.11) 
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where p is the change in length of the resonator 
bring the system back into resonance after the sac 
inserted. <j>. and (j) d arise from the extra phase shi: 
fact that we have a Gaussian beam in the resonat 
plane wave. The permittivity E'=n 2 is then 
solving Eq. (3.7) for n. 

The tangent of the loss angle 5 is given by 

tan5=Q“ : [ (tA+D) / (tA+2k _i sin2 (kd-<|) d ) ] 
where A=n 2 / [n ; cos 2 nkt-d> t ) +sin 2 (nkt-<j) t ) ] 

Q e is the Q factor of the resonator considering 
loss in sample. In open resonators each contrib 
resonant mode losses can be described in terms c 
Q-factors, Qi, which combine in parallel to give 
resonator Q-factor, Q=l/ (eQ^ ) . Let Q d be the measu 
for the resonator containing the sample and Q L th 
Q factor for the resonator containing an idea 
sample of the same dimension and permittivity 
sample. Q L is given by 

Q l =Q 0 (M+1) (tA+d) / (MA+1) D 

where Q d is the measured Q factor of the empty res 
Q e can be calculated by 


required to 
■pie has been 
:t due to the 
or and not a 
obtained by 


(3.12) 

(3 . 13) 

only energy 
union to the 
f individual 
the overall 
red Q-factor 
e calculated 
1 loss free 
as the real 

(3.14) 

coat or. Thus 


i/Qe=l/C,-l/Q, 


(3.15) 
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The above theory assume that the upper surface cf the sample 
is slightly c cr.vex as shown in Figure 5 (a) to match the phase 
front of the beam in the resonator. This means that for the 
plane parallel sample used in practice the small volume shown 
shaded in Figure 5 (b) is unaccounted for the theory. 

phase front 


Fig. 5. Mismatch at sample surface 
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) 2 /4R 
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(3.16) 
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26 


3 . 3 Experimental arrangement 

Developing a mechanism to effectively feed and monitor 
the energy stored in an open resonator is critical in making 
measurements of metallic surface resistivities and dielectric 
constants. The experimental hole-coupling geometry fed by an 
optically pumped far infrared laser is used as shown in Figure 
6. This experiment was performed by Cook ez al.. The 
theoretical analysis of the data is dependent on the fundamen- 
tal mode; higher mode degeneracies must be avoided as we 
mentioned before. 



The experimental apparatus 


The resonate 

r cavity was excited 

by 

the 

output 

from an 

Apcllo Model 122 

FIR laser at NASA 

Lew i 

s Re 

search 

Center, 

operating with a 

Gaussian TEM CC mode 

an 

337 

GHz . 

In this 


arrangement the incident beam and exit beam are separated so 
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that only energy that leaked out from the cavity through the 
coupling hole was monitored. The 2-mm-diameter coupling hole 
should contribute losses of about 2% to the energy stored in 
the Gaussian field in the cavity 7 . Therefore, Eq. (1.46) 
should be multiplied by the factor 0.98 in the theoretical 
calculation due to this coupling hole loss. In cur reference 
cavity both mirrors are aluminum. The radius of curvature of 
the spherical mirror used in measuring metallic surface 
resistivities is 135.9 mrrf. The resonator lengtn at the first 
resonance is just smaller than 135.9 mm (See Table 1) . These 
parameters were calculated by using Eq. (1.10) to limit the 
beam waist w a on the planar mirror to approximately 1.5 mm. 
When measuring properties of dielectric materials, the radius 
of curvature of the spherical mirror was changed to 113.82 mm, 
which decrease the beam waist w Q to less than 1.5 mm. In both 
experimental arrangements, the dimension of the spnencar 
mirror is a little larger than 100 mm, about 2 times the size 
of the Gaussian field on the mirror surface. This large 
curved mirror insures that leakage around the mirror edge can 
be neglected. The planar mirror is masked by a 1/4" (6.4 mm) 
diameter aperture to prevent spurious modes in the cavity. 
This size fits the optimal range given in section 2.2. 
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3.4 Empirical determination of the cavity Q factor 

A measure of the sharpness of response of one cavity to 
external excitation is the Q-factor of the cavity, defined by 
Eq. (1.17) . By conservation of energy the power dissipated in 
ohmic losses is the negative of the time rate c^ change Oi 
stored energy. An initial amount of energy decays away 
exponentially with a decay constant inversely proportional to 
Q and the frequency separation Af between half-power points 
determines the Q of cavity is 

Q=f / Af (3.17) 


Empirically, the Q ractcr ca 
is the length separation of 
resonant length D. To show th 
given by Eq. (1.14). It can 


n be found' 4 from D/AD where AD 
half-power width versus cavity 
is, we use the resonant condition 
also be written as 


f mnq = c [qt-i- (1-m+n) cos "V (R 2 ~D) /R 2 ] /2%D (3.18) 

In the hemispherical resonator, D is very close to R 2 so tnat 
V (R 2 -D> /R 2 much less than 1. If let g= [ (R 2 -D) /R 2 ] ' 4 , the 
expression of ccs“~g in series is 


cos -1 g = ( tc / 2 ) -g- (g 5 / 6) - * (tc/ 2) -g 


(3.19) 
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Then, introducing the fundamental mode TEM ooq into the resonant 
condition, we have 

f=c[q+l/2-g/Jl]/2D (3.20) 

=c [q+l/2-JfW (R 2 -D) /R 2 ] /2D (3.21) 


For V (R 2 -D) /R 2 « 1 and q »1, Eq. (3.21) becomes 


f=C (q+1/2) /2D 


(3.22) 


and resonant 


frequencies of the two half power points are 


f. 2 =f±Af/2 = 0 . 5c [ (q+1/2) / (D*AD) ] 


where the frequency separati 


ation Af is given by 


Af =f ,-f , = Q . 5c (q+1/2) AD/ (D 2 -AD 2 /4) (3.24) 


=0 . 5c (q+1/2) AD/D 2 


(3.25) 


Here we drop the term AD 2 / 4 since D 2 /AD 2 was on the order of 
10 13 in our experiments. According to Eq. (3.17), we have 

Q=Af /f=AD/D (3.26) 

Although it is difficult to make accurate measurements of two 
half-power points D : and D 2 , the half-power width AD=D : -D 2 can 


be measured accurate 
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3.5 Measurement of surface resistivities of metallic samples 

Let Q o1 be the first measured value with mirrors made from 
same metallic material. This serves as a reference to compare 
other sample losses. An attempt was made to give these two 
mirrors similar surface treatments so that the resistivities 
would be the same. Assuming that M=1 in this case, we can 
calculate R cl from Eg. (1.46). Using another different 
metallic sample serving as a planar mirror, we men get Q C2 . 
From Eq. (1.46), it is seen that the ratio of Q C1 to Q 02 is 

Qoi/Q: 2=( m+1 > /2 (3.27) 

which crives us a value of M. Since Q 31 and Q G2 car. oe measures 
in both cases, we therefore use M to calculate the sample 
resistivity R c2 . 

Table 1 shows some sample surface resistivities that were 
calculated from Q values measured at 337 GHz. Note that these 
relative values are based on the literature values from 
aluminum. The literature values are calculated from scaled dc 
values for resistivity 1/So, where 5 and <7 are me classical 
ac skin depth and dc conductivity, respectively. Therefore, 
the literature values should be considered only as rough or 
approximate values. On the other hand, the experiment can be 
considerably improved to limit experimental error on the orcer 
of 5% or lower (see Ref. 1) . 
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Table 1. Surface resistivity of selected metals 


Material 

Q a 

D 

R (measured) 

R (Ref . ) 


10 3 

(mm) 

(10~ 3 O.) 

(10 -3 Q) 

Aluminum 

110 

135.5 

— 

18 8 b 

Copper 

140 

135.0 

107 

14 9 b 

Tantalum 

85 

135.5 

299 


304 Stainless 

30 

133.7 

1191 

i 

97 9 C j 

. 5um thick Au 

96 

135.5 

243 

187 c 


a. From J. D. Cook et al. Rev. Sci. instrum. 62 , 2480 (1991) 

b. From K. J. Button, "Infrared and Millimeter Waves", 
Academic, 1979, Vol.l, p. 237 

c. Calculated from ac conductivity values in R. C. Weast and 
M. J. Astle, "Handbook of Chemistry and Physics", 

CRC, Cleveland, OH, 1975-76, 56th Ed., p. D-171 
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3 . 6 Dielectric permittivity measurements 

The details of the apparatus which operates at 337 GHz 
have been shown in Figure 4 . The sample was placed on the 
planar mirror. The resonant length D for an empty resonator 
can be calculated from a knowledge of the frequency and the 
radius of curvature of the concave mirror. This measurement 
together with the sample thickness t gives us the data 
required for the solving' of Eq. (3.7) for n and hence £' as 
shown in previous section. 

3.6.1 Computational procedure 

Experimental data supplied by Cook et a 1 . are shown m 
Table 2. We use data of crystal quartz, sapphire, and diamond 
no illustrate our technique. 


Table 2. Sample thickness and reduction in resonant length 


Sample 

t ( mm ) 

p (mm) 

Sapphire 

0.508 

0 . 352 

Quartz 

0.659 

0 . 338 

Diamond 

0.250 

0 . 353 
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The sample index n of refraction can be calculated with 
the given data t and p. Note that the correction of the 
sample thickness should be made in calculations according to 
Eq. (3.16). To find the numerical solution of Eq. (3.7) we 
use Math CAD's "solve block" in our calculation. Math CAD 
uses an iterative method to solve equations and inequalities 
in a "solve block". This iterative method starts with guess 
values for the variables to be solved and the constraints and 
it ends with variable values that satisfy the constraints. 
Our "solve block" is listed in Appendix. Table 3 shows the 
results obtained for quartz, sapphire, and diamond. 

3.6.2 Analysis of results 

Equation (3.7) gives multiple solutions or the sample 
index. The "true index" is found among these values by 
comparing these values to those from classical dispersion 
theory. The theory 1 ' applied to a set of lattice oscillators 
relates the dielectric constant £ to the characteristic 
oscillator parameters. In a spectral region where f<< f., for 
all j, the expression for the real part of the dielectric 
constant given by Eq. (3.4) and Eq. (3.5) can be expanded to 
give 
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Table 3. Calculated results for quartz, sapphire and diamond 


Sample 

t (mm) 

p (mm) 

n (guess) 

^val 

t- * ! 

C val i 

t 

Quartz 

0. 659 

0.338 

1 . 5 

1.546 

0.665 ! 

i 

2 . 0 

2.250 

0.664 j 

j 

2 . 5 

2.940 

0.664 | 

3 . 0 

2 . 940 

0.664 | 

3.5 

3 . 628 

0.663 | 

j 

j 

Sapphire 

0.508 

0.352 

2.0 

1 . 635 

0.513 j 
j 

2.5 

2 . 458 

0.513 j 

lj 

3 . 0 

2 . 458 

0.513 || 

ji 

3.5 

3.288 

0.513 lj 

4 . 0 

4 . 125 

0.513 j 

Diamond 

0.250 

0.353 

1 . 5 

— 

i 

2 . 0 

2.53 

0.253 j 
i 

2.5 

2.53 

0.253 'j 
ji 

.'j 

3.0 

— 

J 

ji 

j 

3.5 

— 

•1 

;! 

lj 


a. t val is determined by t+w (t) 2 /4R (t) 

















35 



£ ' =n 2 -K 2 = S c + (S 0 -£..) 

(3.28) 


£" = 2nK= (S c -eJ y 0 X 0 f 

(3.29) 

here 

S 0 =£«+Ij5,=n^ 

(3.30) 


^(So-eJ^ZjSja-Yj) f: 2 

(3.31) 

and 

y 0 =fo 1 (S G -e„)- 1 I j s j y j f; 1 

(3.32) 


where S, is the jth oscillator, f. is the characteristic 
frequency, y. is the dimensionless, frequency-independent 
damping constant, £„ is the high-frequency dielectric con- 
stant, and the S 0 is the zero-frequency or static dielectric 
constant. For transparent substances, such that n 2 >>k 2 , the 
formula for index of refraction derived from Eq. (3.28) is 

n-n 0 + (n;-e.)n'Xf 2 /2 (3.33) 

By plotting experimental data for n vs . f 2 one should obtain 
a straight line according to Eq. (3.33) . The values of n 0 and 
A. 0 may be determined from the intercept and slope of this 
line, provided that n„ is estimated from independent data in 
visible light. Roberts and Coon 16 gave a theoretical discus- 
sion and experiment results at a wavelength range from 4 0(1 to 
300)1 using a grating monochromator for the determination of 
the indexes of quartz and sapphire. Their results, extrapola- 
ted t o the low frequency of about 30 GKz, are: n x =2.102, 
n,=2.153 in quartz and n x =3.060, n,=3.402 in sapphire. Since 
quartz and sapphire are anisotropic materials, we use n x and 
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n, to represent the crystal having its optical axis perpendic- 
ular or parallel to the surface. For low loss materials we do 
not expect any significant dispersion between the low-frequen- 
cy values and the measurements carried out at 337 GHz. Hence, 
we compared our results and the values obtained from the 
literature. Taking these values as reference, we find 2.250 
for quartz and 3.288 for sapphire from Table 3 are the correct 
values of these indexes. The index of diamond as found from 
Table 3 is unique. Van Camp et al. 17 provided the 
literature value n=2.36 for diamond index at the optical 
wavelength near 4|i. 


Table 4. Comparison of refractive indexes 


Material 

n val (Cal.) 

n (Ref.) 

Variation 

Quartz 

2 . 25 

2 . 102 

7% 

2.153 

5% 

Sapphire 

3.28 

3.060 

7% 

3.402 

4% 

Diamond 

[| 

2 . 53 

2.36 

7% 


Table 4 compares our results with the reference values. 
The calculated n val values for quartz and sapphire are based on 
the samples with unknown optical axis and the corresponding 
reference values are between n x and n, . However, Table 4 shows 


However 
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a consistency with literature values cf about 7%. 
open resonator and Math CAD computational proce 
useful method for determining the permittivities c 
samples. It is important to point out that t 
modified in anisotropic samples. Emert 18 and Bha 
shown that the maximum effect of dielectric anis: 
estimated by substituting t" for t' in Eqs . (3. 

where 


t" = (n 2 /nj) t' = (n,,/nl) [t+w (t) 2 /4R(t) ] 

With this change, more accurate measurements c 
indexes for anisotropic samples may be obtained. 


In general, 
cures give a 
f dielectric 
should be 
walkar 19 have 
:cropy can be 
") to (3.11) 

(3.34) 

: refractive 
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3.7 Measurement of loss tangent in some transparent 
substances 

Loss-tangent measurements are basically dependent on the 
measured Q values in the resonator. It is desirable uhat Q<j 
should be, within an order of magnitude, comparable with Q 0 in 
low loss materials. If it is not, one should be able to carry- 
out very precise and repeatable Q measurements. 

3.7.1 Results 

Empirical Q measurements determine D/AD by defecting the 
half-power points . Our calculated results are based on the 
data from the experiment performed by Cook et al . at u37 GHz. 
The samples are the same as we used in dielectric permittivity 
measurements. Figure 7 shows comparative cavity scans for some 
dielectric samples. The measured values of half-power width 
are given in Table 5. 



Power, Arb. Units 


39 





J;r^ \}mhs 


'^jdi, /'tiW'LMn jKtAt^ J! X ( O ) 


100 


30 


60 


40 


20 


Change in mirror separation, [i 

Fig. 7. Sample spectra taken with a scanning hemispherical 
open resonator for empirical cavity Q measurements 
(a) empty (b) loaded with quartz 

(c) loaded with diamond (d) loaded with sapphire 



Table 5. Q values for the empty and sample-loaded resonator 


Sample 

D (mm) 

AD/ ((I) 

AD. 3 (|1) 

Q 0 C ( 10- 1 ’ 

Q« a <i0 3 ) 

Quartz 

113 .707 

1.25 

1.50 

mm 

7 6 

Diamond 

113.707 

1.25 

2.25 

■i 

j 


113 .707 

1.25 

1.56 

mm 

65 ! 


a. Half-power width for the empty resonator 

b. Half-power width for the loaded resonator 

c. Q-factor for the empty resonator 

d. Q-factor for the loaded resonator 

Using Eqs . (3.12), (3.14) and (3.15), we calculated the 
less tangents shown in Table 6 . 


Table 6. Loss tangents of Quartz, Diamond and Sapphire 


Sample 

t (mm) 

p (mm) 

n ref erer.ee 

Q.CLC 3 ) 

BB 

s^uart z 

0.659 

0.333 

2 . 153 


0.035 | 

j 

Diamond 

0.250 

0 .353 

2.36 

1347 

0.154 ■ 

■Httfl 

0.508 

0.352 

3 . 402 

230 

0.654 

! 


The resonator cavity is extremely sensitive to small 
misalignments in either mirror. Therefore, it is critical that 
the translational stage be firmly mounted in such a manner 
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that forward and backward movement will not misalign the 
cavity over multiple runs. This misalignment presented itself 
in this experiment because the intervals between the scaled 
mirror separations shown in Figure 7 are not linear at certain 
rate of cavity scan. It is estimated that errors in determin- 
ation of AD were about 10% so that the uncertainty in tan5 
was found to be 25% due to this method or measuring the Q— 
f actor and other losses in the cavity. 

3.7.2 Discussion 

Mirrors used in the resonator are made from aluminum. 
Using the measured values of surface resistivity of aluminum 
found in Table 1, we can obtain a theoretical empty Q-factor 
from Eg. (1.46). This Q value is about 9.2xI0 4 in the 
resonant length of 113.707 mm. The empty Q-factor Q 0 measured 
t>y the experiment and shown in Table 5 is uhen found to agree 
with the theoretical estimates. We expect to maxe measure- 
ments of Q-factors of magnitudes of about 10 2 to better than 
10% . 

Except for the uncertainty of this measuring method, 
another factor affecting the measurement of loss tangent is 
the sample thickness. It has been shown by Jones (see Ref. 
13) that possible errors are minimum when the sample thickness 
is an intecral number of a half wavelength ihick. Uni ortunat 
ely, none of our samples has a thickness like this because of 
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limited availability. It is expected that the loss tangent of 
sapphire has the least error since its thickness is closer to 
a half wavelength than the other two samples. 

Although the measurement was somewhat inaccurate due to 
the problem mentioned above, the use of the resonator in the 
submillimeter-wave region is a feasible method to measure loss 
tangent. This is the first demonstrated use of a far-infrared 
laser as a dependable research tool in this context. 



4 . Summary 


In this work we nave extended the derivation of theoreti- 
cal formulas for Q-factors of open resonators by using scalar 
theory and Gaussian beam model for small sample. This work is 
used in determining experimental Q-factors ottained from 
measurements of metallic surface resistivities anc dielectric 
constants . It is important because these measurements are 
dependent on the accurate measurement of the cavity Q. To 
avoid spurious oscillating modes in the resonator, use of a 
small sample is essential. The modified theory, dealing with 
the small planar mirror of the open resonator, indicates thcu 
the existing theories can still be used if the sample radius 
is larger than 1.5w 0 . 

This work justifies the use of small sample arrangements. 
Using hole coupling techniques in the open resonator at hign 
frequencies of about 337 GHz with FIR lasers is a new approach 
in measuring metallic surface resistivities anc dieiectnc 
constants. Although the response time and sensitivity or the 
experimental system needs to be improved, it extends applica- 
tions of the open resonator into the submillimeter-waveiengtn 
region . 
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